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ABSTRACT 

Using high-resolution simulations within the Cold and Warm Dark Matter models 
we study the evolution of small scale structure in the Local Volume, a sphere of 8 Mpc 
radius around the Local Group. We compare the observed spectrum of mini-voids in 
the Local Volume with the spectrum of mini-voids determined from the simulations. 
We show that the AWDM model can easily explain both the observed spectrum of 
mini- voids and the presence of low-mass galaxies observed in the Local Volume, pro- 
vided that all haloes with circular velocities greater than 20 km/s host galaxies. On 
the contrary within the ACDM model the distribution of the simulated mini- voids re- 
flects the observed one if haloes with maximal circular velocities larger than 35 km s~^ 
host galaxies. This assumption is in contradiction with observations of galaxies with 
circular velocities as low as 20 km/s in our Local Universe. A potential problem of 
the AWDM model could be the late formation of the haloes in which the gas can be 
efficiently photo-evaporated. Thus star formation is suppressed and low-mass haloes 
might not host any galaxy at all. 

Key words: cosmology: large-scale structure of Universe, voids, dark matter, wdm 
paradigm ; galaxies: luminosity function, rotational velocities 



1 INTRODUCTION 

At present the standard model of cosmology is a fiat Fried- 
mann universe whose mass-energy content is dominated by 
a cosmological constant (the A term), a Cold Dark Mat- 
ter (CDM) component and baryons. This ACDM model is 
characterized by only a few parameters which are deter- 
mined already with high precision. The ACDM model de- 
scribes structure formation at large scales very well, how- 
ever it fails on small scales: the standard model predicts 
much more small scale structure than observed. This con- 
cerns the number o f dark matter sub-h a loes inside galactic - 
sized host haloes (|Klvpin et al.l 1 19991 : iMoore et all ligogl ) 
as well as th e number of d warfs in low density regions of 
the universe (|Peeblesil200lh . A better understanding of the 
physics of structure formation on small scales, in particu- 
lar of the correc t modeling of baryo nic physics, coul d solve 
these problems l|Hoeft et al.,.2006. : .Crain et al.ll2007l ). How- 
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ever, any non-baryonic physics that reduces power on small 
scales compared to the standard model will also improve the 
situation. It is well known that warm dark matter acts in 
this direction by erasing power at short scales due to free 
streaming. 

Recently iTikhonov fc KlvpinI (|2009l ) studied the spec- 
trum of mini-voids in the distribution of galaxies down to 
Mb = —12 in a sphere of radius 8 Mpc around the Milky 
Way (the Local Volume - LV). They compared the spectrum 
of mini-voids obtained in LV-candi dates selected in nu meri- 
cal simulations based on WMAPl (jSpergel et al.ll2003l ) and 
WMAP3 l|Spergel et al.l bOOTl ) cosmologies and concluded 
that it matches the observed one only if haloes with circular 
velocities 14 > 35 — 40 km /s define the mini-voids. Thus 
iTikhonov fc KlvpinI ||2009D conclude that ACDM haloes 
with Vc < 35 — 40 km/s should not host galaxies brighter 
than Mb = —12, in fairly agreement w ith theoretical ex- 
pectations (|Hoeft et al.ll2006l : lLoebll2008t ). In the LV, how- 
ever, a substantial number of quite isolated galaxies has been 
observed with a total of reg ular (circular) and chaotic mo- 
tions well below 35 km s"^ (|Begum et al.ll2008t ). These ob- 
served galaxies point towards the same overabundance prob- 
lem that the ACDM model has in case of Milky Way satel- 
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lites. If haloes with circular velocities Vc < 35km/s would 
host galaxies their total number is predicted to be too large 
compared to data. 

In Warm Dark Matter (WDM) cosmologies the num- 
ber of low-mass haloes drops down when their mass is be- 
low a critical damping scale wh ich depends on the mas s 
of the WDM particles ( see e g lAvila-Reese et al.l l)200lf ). 
iBode. Ostriker. fc Turok (200l|), lKnebe et al.l (|2002h. In this 
work, we will use high-resolution N-body simulations of dif- 
ferent dark matter candidates, namely CDM and WDM, 
in order to study the properties of LVs in these models 
and compared with the real LV observations. In particu- 
lar, we will compare the void functions in simulated "LV- 
candidates" with observations in order to derive possible 
conclusions about the nature of dark matter particles. 

Unfortunately, little is known about the nature of the 
dark matter. One of the promising dark matter candidates is 
the gravitino for which three possible scenarios are imagin- 
able: Warm (super WIMP) or mix ed cold/warm or cold grav- 
itino dark matter (jSteffenl |2006| ). Another possible candi- 
date are sterile neutrinos which exhibit a significant primor- 
dial velocity distribution and thus damp primordial inhomo - 
geneities on small scales (|Abazaiian fc Koushiappa3l2006l ). 
Experiments are far away from the scales which could probe 
directly the properties of the particles of which the dark 
matter is made of, so that the free streaming length can- 
not be predicted. However, the free streaming length can be 
associated with the particles mass for which astrophysical 
constraints can be determined. 

There is a number of observational constraints on the 
mass of the WDM particles. Using HIRES dat a, a lower 
limit of mwDM = 1.2keV has been reported by IViel et al.l 
which increases to 4.0 keV when a combination of 
SDSS and HIRES data was used. From the power spec- 
trum of SPSS galaxies, a lowe r limit of 0.11 keV has been 
reported by [AbazaiianI (|2006l ) which also increases to 1.7 
keV if the SDSS Lya forest is considered and further in- 
creases to 3.0 keV if high-resolution Lya forest data are 
taken into account. An independent estimation of lower lim- 
its for the mass of the WDM particle based on QSO lens- 
ing observations are in a greement w ith the former results 
( Miranda fc Macciol [2OO7I ') . Recently. Iviel. Colberg. fc KimI 
( 20081 ) have shown that the properties of the void popula- 
tion seen in the Lymana forest spectra depend on the under- 
lying power spectrum and, therefore, these voids may also 
constrain the properties of the dark matter. We have run a 
suite of WDM simulations with particle masses ranging from 
3 to 1 KeV. However, in this paper we will consider only the 
results from simulations done with the lighter WDM par- 
ticles (1 KeV). Results from the other simulations will be 
reported elsewhere. 

In this paper we will focus on the comparison between 
the observed spectrum of mini-voids in the Local Universe 
with the simulated ones in ACDM and AWDM models. The 
structure of this paper is as follows, in Section 2 we describe 
briefly our numerical simulations, the halo mass functions 
and discuss the reliability of the subsequent analysis. In Sec- 
tion 3 we describe the observational data sample used here. 
In Section 4 we discuss our void finding algorithm. In Section 
5 we present results from AWDM and ACDM simulations 
in terms of spatial distribution and dynamics. We conclude 
this paper with a discussion of the results found. 



2 SIMULATIONS 

For our numerical experiments, we have assumed a spatially 
fiat cosmological model which is com patible with the 3rd 
year WMAP data (|Spergel et al.ll200'i1 ) with the parameters 
h = 0.73, n„ = 0.24, Qiar = 0.042, = 0.76, the nor- 
malization (Tg = 0.75 and the slope n = 0.95 of the power 
spectrum. A computational box of 64ft.~^Mpc on a size was 
used for all simulations. For ACDM a random realization of 
the initial power spectrum computed from a Boltzmann code 
by W. Hu, was done using 4096'^ N-body particles. After ap- 
plying the Zeldovich approximation to this set of particles, 
we reduced the particles to a total number of 1024^, which 
corresponds to a mass per particle of 1.6 x 10''h~^MQ. We 
used the original displacement field on the 4096^ mesh to 
derive the initial conditions (i.e positions and velocities) of 
the final 1024'^ particles. In this way, we can increase the 
mass resolution of all, or a particular area of the simulation 
up to 64 times more using zooming techniques. This is spe- 
cially useful to study the problem of fake haloes in WDM 
(see below). The same random phases have been used in 
both simulation, so we ensure that the same structures will 
be formed in both models. The initial redshift of the simu- 
lations was set to 50 in all cases. 

To generate initial conditions for WDM, we used a 
rescaled version of the CDM power spectrum using a fit- 
ting function that approximates the transfer function asso- 
ciated to the free stre aming effect of WDM particles with 
mwDM = 1 — 3 keV l|Viel et al.ll2005l ). The power spectra 
of the three models are shown in Fig. [T] 

The effects of the velocity dispersion of the WDM par- 
ticles have been neglected in these simulations. They would 
introduce a certain level of white noise in the initial condi- 
tions t hat would contribute to the growth of th e small scale 
power (|Colm. Valenzuela. fc Avila-Reesell2008l ) that lead to 
the formation of spurious halos. At the same time, the rms 
value of the random velocity component for 1 kev parti- 
cles at the starting redshift of the simulations (z=50) is of 
the order of 2 km/s which is much smaller than the typical 
velocities induced by the gravitational collaps e of the mini- 
mum resolved structures in our simulation (see lZavala et al.] 
(,2009 ) for more details). 

These simulations were run with the purpose in mind of 
studying the gross structures of the nearby Universe, such 
as the Local Supercluster (LSC) and the Virgo cluster, i.e, 
to use them to do Near Field Cosmology. This could be 
achieved by imposing linear constrains coming from obser- 
vations to the otherwise G aussian random fie l d of d ensity 
fiuctuations. We used the iHoffman fc RibakI (|l99ll ) algo- 
rithm for imposing the constrains on large scales. A more 
detailed description of the simulations and their use for do- 
ing Ne ar Field Cosmolog y will be reported in a forthcoming 
paper (| Yepes et al.ll2009l ). However, the constrained nature 
of the simulations is ignored here and they are treated as 
standard random simulations. 

The simulations have been per formed using t he TreePM 
parallel N-body code GADGET2 (IS pring ell lioOsI ). The spa- 
tial resolution (Plummer softening length) was set to the 
maximum between 1.6/i~^kpc comoving and 0.8/i"^kpc 
physical at all redshifts. 

In an N-body simulation, the computational box size 
(kp) and the number of particles (kMy) determine the fil- 
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Figure 1. The dimensionless theoretical power spectra for the 
models considered here: CDM (black), and WDM with particle 
masses: 3 keV (blue), and 1 keV (red). The dotted lines corre- 
spond to the fundamental mode (box size 64/i~^Mpc) and the 
Nyquist frequencies. 



Figure 2. Cumulative EOF halo mass functions with more than 
20 particles for the ACDM (dashed line) and AWDM (solid 
line) models. From left to right the mass functions are shown 
at redshifts 2: = 8, 6, 5, and 0. The dotted line marks Miim = 
3 X IO^/i-^Mq. 



tering scales of the theoretical power spectrum. In Fig. [T] 
we represent these two scales for our simulations. As can be 
seen, the Nyquist frequency is well below the exponential 
cut-off in the power spectrum for WDM with raw dm = 3 
keV. Therefore, we cannot expect too much difference be- 
tween the halo mass function of ACDM and this particu- 
lar AWDM model, since both are dominated mainly by the 
lack of mass resolution. Therefore, we decided not to run 
the 3 keV WDM simulation at this resolution. Rather, we 
used it to study the excess of satellites in Milky Way type 
haloes, in which re-simulations of haloes with full resolu- 
tion (i.e. 409 6^ particles) have been used for the 3 models 
ijYepes et al ] r2009 ). Here we concentrate on the simulation 
with 1024^ WDM-particles of 1 keV mass. 

In order to identify haloes in the simul ation we have 
used the AMIGA Halo Finder (AHF) (Kno Umann fc Knebd 
l2008l ). The halo positions are first identified from the local 
maxima of the density field that is computed from the par- 
ticle distribution by Adaptive Mesh Refinement techniques. 
The AHF identifies haloes as well as sub-haloes. It provides 
the virial masses, virial radii and circular velocities of iso- 
lated haloes and some related quantities for sub-haloes. We 
use the maximum circular velocity 14 to characterize the 
haloes because the circular velocity can be easier related 
to observations than the virial mass. For reference, haloes 
with Vc ~ 50 km s~^ have a virial mass of about 10^'^ 
and haloes with 14 = 20 km s^^ have a virial mass about 
IO^Mq. We re-scaled all data (coordinates and masses of 
haloes) to "real" units assuming jJp = 72 km/s/M pc, which 
is the value assumed bv lKarachentsev et al.l l|2004l ) in the ob- 



served LV and which is almost identical with value assumed 
in the simulations (_ffo = 73 km/s/Mpc) 

It is known that in WDM simulations small mass haloes 
can arise from the artificial gravitational growth of Poisso- 
nian fiuctuati ons that lead to an un physical fragmentations 
of filaments (|Wang fc White! |2007| ). These "fake" haloes 
show up typically as tiny haloes regularly aligned along fila- 
ments. In Figure[T]one can clearly see the exponential cut-off 
of power in the AWDM models. For the IkeV WDM model 
the maximum power is reached at fcpoak ~ 3.7/iMpc~^. 
Above this wave number the linear power spectrum contains 
very little power so that discreteness noise influences the 
growth of small-scale non-linear structures. Below a charac- 
teristic mass Afiim, which depends on the mass resolution 
in the simulation, a spurious fragmentation of filaments oc- 
curs which leads to an upturn of the halo mass function 
at small scales (Fig. [gjl. Following IWang fc White! (|2007l ) 
we have estimated for our simulations this limiting mass as 
Afihn « 3 X lO^/i"^M0. 

Figure [2] shows the evolution of the cumulative mass 
functions of halos found by the friend-of-friends algorithm 
in the ACDM and AWDM simulations. The upturn of the 
AWDM mass function can be seen at any redshift between 
2 = 8 and 2 = and it is always at the same position slightly 
left of the limiting mass. The dotted line marks the limiting 
mass 3 x lO^/i~^M0 for 2 = 0. This points to a complex 
process of the formation of spurious fake-haloes. It seems 
that they do not form only at at certain epoch but instead 
appear all the time after the formation of the first haloes. A 
more detailed numerical study of this problem is definitively 
needed, but it is outside of the scope of the present paper. 
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Since we see an upturn of the mass function only at 
~ 1 X 10®/i~^Mq we consider Mum as a rather conservative 
value. In terms of circular velocity Muptum corresponds to 
haloes with Vmax 20km/s which are candidates to be spu- 
rious. The prime objective of the paper is the study of the 
statistical properties of voids, and it is therefore important 
to establish that these are not affected by the spurious halos. 
Figure [5] below proves that this is indeed the case. 



3 OBSERVATIONAL DATA: LOCAL VOLUME 

Over the last decade Hubble Space Telescope observations 
have provided distances to many nearby galaxies which 
where measured with the tip of the Red Giant Branch 
(TRGB) stars. Special sea rches for new nearby dwarf galax- 
ies have been undertaken iKarachentsev et al.l (|2004 ). Also 
speci al cross-check has been don e with IR and blind HI sur- 
veys l|Karachentsev et al.l l2007t ) . The distances were mea- 
sured from TRGB stars, Cepheids, the TuUy-Fisher rela- 
tion, and some other secondary distance indicators. Since 
the distances of galaxies in the Local Volume are measured 
independently of redshifts, we know their true 3D spatial dis- 
tribution and their radial velocities. At present, the sample 
contains about 550 galaxies with distances less than 10 Mpc. 
A substantial fraction of the distances have been measured 
with accuracies as good as 8-10% ( Karachentscv et al. 2004). 

Within the Local Volume dwarf systems have been ob- 
served down to extremely low luminosities. This gives us 
a unique chance to detect voids, which may b e emp ty of 
rather faint ga laxies. Tikho nov fc KarachentsevI (|2006l ) and 
iTikhonov fc K lvpin (2009) analyzed ne arby voids. Here we 
use the same observational data as in ITikhonov fc KlvpinI 
2009[) - the updated Catalog of Neighboring Gala^des 



Karachentscv et al.ll200'j . private communication). We an- 



alyze a volume-limited sample of galaxies with absolute 
magnitudes Mb < —12 within 8 Mpc radius. The over- 
all completeness o f the sample has been discussed in 
ITikhonov fc KlvpinI (|2009l ). 



3.1 Circular vs. rotational velocity 

In order to compare results from an observational galaxy 
sample with dissipationless N-body simulations one has to 
assume certain hypotheses about how galactic properties are 
related with their host dark matter halos and, in particular, 
how the observed rotational velocity of a galaxy is related 
to the circular velocity of its dark matter halo. 

A number of different techniques has been developed to 
associate galaxies with simulated DM haloes. The most well 
known te chniques are the semi-analytical models (see the 
review of iBaughj l|2006l l). An alternative approach is based 
on th e conditional luminosity function (jvan den Bos ch et al.l 
|2004 ). Based on the merging history of the DM haloes and a 
set of physically motivated recipes the semi-analytical mod- 
els place galaxies into the DM haloes. A number of free pa- 
rameters of these recipes are calibrated by the observation 
of galaxies. A dark matter halo may host a luminous galaxy 
and a number of satelittes. Although modern high resolu- 
tion DM simulations resolve sub-haloes most of the semi- 
analytical models do not follow haloes which are accreted 
by larger ones but rather model the transformation of these 



accreted satellite s along their orbits. R ecently, lMaccio' et all 
Hooi) (see also IKoposov etai] Hooi)) claimed that semi- 
analytical models can solve the lon g standing problem of 
missing satellites (|Klvpin et al.|[l999l ) within the ACDM sce- 
nario. The basic idea behind this solution of the problem is 
that the haloes which host a galaxy of a given measured 
rotational velocity are more massive than expected by the 
direct association of rotational velocity and the haloes max- 
imum circular velocity. Since more massive haloes are less 
frequent the problem of missing satellites is solved. Thus 
it is solved at the expense of a predicted large population 
of dark satellites the observation of which remains another 
ope n problem in cosmology. 

IStoehr et all (|2002l ) were one of the first who argued that 
due to gravitational tides the total mass associated with 
dwarf spheroidals might be much larger than assumed. On 
the contrary, Kazantzidis ct al. (2004) argue that tidal inter- 
action cannot provide the mechanism to associate the dwarf 
spheroidals of the Milky Way with the mas sive substructures 
predic ted by the ACDM model. Recently, iPefiarrubia et ah! 
|200») analyzed the dynamics of the stellar systems in local 
group dwarf spheroidals and came to the conclusion that the 
associated dark matter circular velocity is about three times 
higher than the measured velocity dispersion of the stars. 

All the previous arguments are mainly based on the as- 
sumption that the small galaxies under consideration are 
satellites of a more massive galaxy. However, in the analy- 
sis we report in this paper we will consider isolated dwarf 
galaxies for which these arguments are not necessarily valid. 

In what follows we argue in favor of a direct association 
of the observed rotational velocity of a isolated dwarf galaxy 
to the maximum circular velocity of a simulated low mass 
DM halo. Our arguments are based on the comparison of 
the beam size of the observations with theoretical profiles of 
the DM haloes. 

The IKarachentsev et al.l l|2004l) galaxy catalog contains 
for many of its galaxies the 21 cm HI line width at the 50% 
level from the maximum (W50) as well as the apparent ax- 
ial ratio {b/a). Using these informations ITikhonov fc KlvpirJ 
(^2009) estimated the rotational velocity of the galaxy as 
Kot = W^5o/2Vl - (&/«)'• They argue that a large fraction 
of the HI line width is likely produced by random motion of 
the order of 10 km s~^ but did not subtract those random 
motions. Their conclusions were based on a comparison of 
this rotational velocity with the circular velocities of the 
haloes. We follow their suppositions and associate the ro- 
tational velocity K-ot with the peak circular velocity of the 
dark matter halo. In order to give additional arguments in 
support of this assumption we compare now circular veloc- 
ity profiles of dark matter haloes with the HI radius within 
which the observations are performed. 

In Figure |4] we plot the NFW velocity profiles up to 
4_Rs for 7 haloes with masses between 5 x lO**/i~^M0 and 
2 X lO^°/i~^M0. Here, we assu med concentratio ns following 
the fitting form ula provided by Netq et aP (|2007l ) which is in 
agreement with iMaccio et a"iTl|2007l ). Now we want to com- 
pare these circular velocities with the measured rotational 
vel ocities of near b y dwa rfs. We use the FIGGS galaxy sur- 
vey |BegumelaL| (120081 ) with Mb > -14.5 which is located 
inside the Local Volume. In this survey, measurements of 
W50 together with the HI diameters at a column density 
of « 10^^ atoms cm~'^ and the inclinations of optical images 
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Figure 3. The AHF cumulative halo mass functions at redshift 2 = (left) and the corresponding velocity function (right) for the 
ACDM (solid line) and AWDM (dashed-dotted line ) simulations. The two dotted lines indicate the behavior of the mass function 
(velocity function) for haloes which are very likely fakes and those which should be real ones. We mark the intersection point as the mass 
limit below which fake and real haloes are mix ed in the WDM simulation. For comparison, we also indicate the value for this mass limit 
calculated according to IWang fc Whitg l|20o3). 
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Figure 4. NFW circular velocity profiles for DM haloes with 
masses between 5 X IO^H-'^Mq and 2 X 10^°h-'^MQ. Red stars: 
observed Vrot vs HI diameters taken from the FIGGS LV galaxy 
sample. 



(iopt) are provided. We estimated the rotational velocities as 
Vrot = VK5o/(2 X sin(jopt)) for galaxies with iopt > 40° and 
plotted them vs. the HI radii on top of the NFW profiles in 
Figure O 

One can clearly see that most of the observational 
points are located on the flattening parts of the NFW veloc- 
ity profiles with circular velocities below 30 km s~^. Since 
dwarf galaxies are expected to be dark matter dominated, 
the association of the observed rotational velocities at the 
HI radius with the maximum circular velocity of DM haloes 
is in our opinion a valid approximation. 



4 VOID FINDER 

Clusters, filaments and voids constitute the main compo- 
nents of the large-scale structure of the Universe. Soon after 
the discovery of the first cosmic void in the Bootes con- 
stellation (Kirslmc r et. al,, 19811 ) it became clear that voids 
in galaxy distribu tion are the natural outcome of gravita - 
tional instability (|Peeblesl [l98^ : iHoffman fc ShahamI [l983 ) 
During the last two decades the void phenomenon has been 
discussed extensively both from a theoretical as w ell as a 
observational point of view (see e.g. jPeeblesI (1200 ll ). There 
are many different possibilities to detect voids. A thorough 
review and comparison of differe nt algorithms of void find- 
ers has been recently presented in lColberg et al.l l)2008h . The 
void finders can be classified according to the method em- 
ployed. Most are based on the point distribution of galaxies 
or DM haloes and some on the smoothed density or potential 
fields. Some of the finders are based on spherical filters while 
others assume no inherent symmetry a n consider voids a cru- 
cial co mponent of the cosmic web (e.g. iForero-Romero et al.l 

ifiS)). 

Our void finder detects empty regions in a point dis- 
tribution. These regions are not necessarily spherical. Our 
method is not sensitive to the total number of objects since 
the objects are mostly located in high density regions. But 
it is very sensitive to the appearance of objects in low den- 
sity regions. In fact, any object detected in the middle of a 
void would divide it into two vo ids. A detailed description 
of the algorithm can be found in iTikhonov fc KarachentsevI 
(200^). Our void find er is s omewhat similar to th at used by 
lEl-Ad fc PiranI (jl997l l and ICottlober et all (j2003l ). Namely, 
we place a 3D mesh on the observational or simulated sam- 
ple. On this mesh, the identification of voids starts with the 
largest one. At first, we search for the largest empty sphere 
within the considered volume. The voids are supposed to be 
completely inside the geometric boundaries of the sample. 
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therefore, for a given grid point the radius of the largest 
empty sphere is determined by the minimum of the dis- 
tance to the nearest galaxy (or halo) and the distance to 
the boundary of the sample. 

The consecutive search for empty seed spheres with de- 
creasing radius, Raecd, is combined with their subsequent ex- 
pansion by means of addition of spheres whose centers are lo- 
cated inside the given void and whose radius is _Ri > k-Rsccd 
where we assume k = 0.9. Thus, the spheres added to the 
void intersect with it by more than 30% of the volume. At 
any search step for an empty seed sphere, all previously iden- 
tified voids are taken into account, i.e. void overlapping is 
not allowed. The search for voids continues as long as -Rscod is 
larger than a threshold of 0.5 Mpc. This algorithm is simple, 
flexible and appropriate for our definition of voids as regions 
completely free from galaxies of certain luminosities. 

With our assumption, k — 0.9, the voids adequately 
represent the empty regions between galaxies and at the 
same time preserve a sufficiently regular (elliptical) form. 
For k — 1, voids would be strictly spherical. When k de- 
creases, the voids begin to penetrate into smaller and smaller 
holes, and the shapes of voids become more irregular. As- 
suming a very small k, one single void (the first one) would 
fill most part of the sampled volume. The value of the co- 
efficient k = 0.9 used in this paper is a compromise, chosen 
after empirical detection of voids in distributions of points 
with different properties. The voids turn out to be well sep- 
arated one from another and they can be approximated by 
triaxial ellipsoids. 

Once we have identified the voids, we calculate the cu- 
mulative void function (VF) of their volumes. We applied 
also a spherical variant (fc — 1) of our void finder and found 
that the assumption of sphericity of voids does not affect 
our results. 



5 RESULTS 

5.1 Local Volumes in the simulations 

In the following we use the same criteria as 
iTikhonov fc KlvpinI l|2009l) to select "LV-candidates" 
from our simulations. These are spheres of radius 8 Mpc 
which are characterized by the main features of the ob- 
served Local Volume galaxy sample. In our simulations the 
LV-candidates must be centered on a DM-halo with virial 
mass in the range 1 x 10^^ Mq < rrivir < 3 x IO^^Mq. 
This is our definition of a Local Group candidate. We do 
not require that this local group consists of two haloes 
comparable in mass because the appearance of two haloes 
instead of one halo of the same total mass does not change 
the dynamics of matter outside the Local group. We also 
do not require that the Local group should be at the 
same distance from a Virgo-type cluster as the observed 
one. Instead, we impose an overdensity constraints on the 
distribution of mass in the LV, which we take from the 
observed LV sample of galaxies. Below we summarize the 
conditions, which were used to find Local volumes in our 
simulations: 

(i) We put a sphere of radius 8 Mpc on a Local group 
candidate. 



(ii) No haloes with m„ir > 2 x IO'^^Mq are inside this 
sphere. 

(iii) There are no haloes more massive than 5.0 X 10^^ M0 
in a distance between 1 and 3 Mpc. 

(iv) The centers of the spheres (the Local Group candi- 
dates) are located at distances larger than 5 Mpc one from 
the other. 

(v) The number density of haloes with 14 > 100 km s~^ 
inside this sphere exceeds the mean value in the whole box 
by a factor in the range 1.4 — 1.8. 



We found 14 LV-candidates in the AWDM simulation 
which fulfill these criteria and 9 in the ACDM simulation. 
The reason for the diff'erence is the tight requirement on the 
number density of massive halos. This means that in total 
between 18 and 23 massive halos must be in this sphere. 
Thus small differences between the position of a halo or 
its circular velocities in the ACDM and AWDM simula- 
tions may already lead to a violation of this criterion. Since 
the large scale structure in both simulations is quite sim- 
ilar we decided to use in the ACDM simulation the same 
LV-candidates as in the AWDM simulation. All these LV- 
candidates fulfill the first four criteria, 7 fulfill also the over- 
density criterion, 5 are 10 % above or below (what is equiva- 
lent to 2 additional or missing halos with Vc > 100 km s^^) 
and two miss 5 halos. However, all these five halos are in 
the range 90 km s~^ < Vc < 100 km s~^. The advantage 
of our choice is that we can compare the properties of the 
same regions with the same large scale environment in both 
simulations. 

Halo finders like FOF or AHF can find haloes down 
to 20 particles, however to determine the internal proper- 
ties of haloes one needs much more particles. Moreover, 
the sample is not necessarily complete down to 20 parti- 
cles per halo. There is a large scatter in the Vc vs. Npart 
relation which reduces substantially if one considers only 
isolated haloes. Haloes with Npart ~ 100 have circular ve- 
locities of Vc ~ 20km/s. In our simulation this limit co- 
incides with the upturn of WDM velocity-mass function 
due to fake haloes. As an example for the existence of fake 
haloes in AWDM simulations Figure [5] shows the same LV- 
candidate in both models: top row for AWDM bottom row 
for ACDMIn the top, right plot one can clearly see filaments 
of regularly aligned spurious fake-haloes which extend across 
the voids and decrease their sizes. These are haloes with very 
low circular velocities. In the top left plot we have marked 
haloes with 15 < Vc ^ 20 km s~^ by crosses and halos with 
Vc > 20 km s~^ by filled circles. One can clearly see the dif- 
ference to the right plot: Haloes with Vc ^ 20 km s~^ (filled 
circles in the left plot) do not appear at all in artificial chains. 
Haloes with 15 < Vc ^5 20 km s"'^ only slightly modify the 
overall distribution and do not fill voids. Thus, they do not 
change substantially the void statistics. We did not find any 
such artificial regular filament of tiny haloes when we con- 
sider only halos above our limit of reliability Vc = 20km/s. 
We conclude that our analysis is not likely to be affected 
by spurious haloes as long as we consider halo samples with 
circular velocities Vc > 20 km s~^ and might be slightly af- 
fected for samples with Vc > 15 km s~^. 
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Figure 5. Distribution of haloes in a slice of 4 Mpc thickness through the center of an LV-candidate. The upper panel shows the AWDM 
model, on the left plot haloes with Vc ^ 20km s~^ are marked by filled circles and haloes with 15 < Vc < 20km s~^ by crosses. On the 
right plot all haloes are shown. The lower panel shows the distribution of the corresponding ACDM haloes. 



5.2 Properties of voids 

The Local Volumes contain a large number of mini-voids 
with radii Rv smaller than about 4 Mpc. We characterize 
these mini- voids by the cumulative volume of the mini- voids 
Vvoid{R > Rv) normalized to the Local Volume Vlv- In 
Figure |6] we compare the cumulative volume filling fraction 
(VFF) of mini-voids (Koid(i? > R.)/Vlv) of the 14 LVs 
found in the AWDM simulation with the observed VFF for 
Mb < — 12 galaxies. The figure indicates that the VFF in 
the distribution of haloes with Vc > 20 km s~^ fits well to 



the observed one. In fact, the observed spectrum of local 
mini-voids is reproduced within the AWDM model almost 
over the full range of void volumes. The size of the large 
voids reduces slightly if we assume that haloes with Vc > 
15 km define the voids. With the present mass resolution 
we cannot exclude that this reduction is due to fake haloes. 
In the luminosity range —13.2 < Mb < —11.8 only three 
field gala xies have been observed w ith Vrot < 15 km s"'^ 
(table 3 of iTikhonov fc KlvpinI (|2009| )). So far, the existence 
of such very slowly rotating (low mass) galaxies is not in 
contradiction with the assumption of WDM. 
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The VFF obtained from the ACDM simulation matches 
observations only if we assume that objects with Vc > 
35 km s"^ define the local mini- voids (Figure[7j. The ACDM 
VFF is well above the data, if we assume Vc > 50 km s~^ for 
the void definition and it is far below the observational data 
if we assume Vc > 20 km s"'^, i.e. in this case large mini-voids 
are missing in the numerical realizations of local volumes. 
Thus, our ACDM simulation is in agreement with observa- 
tional data if haloes with Vc < 35 km s~^ do not host any 
galaxy. This statement is independent of the volume under 
consideration. As a test of reliability we repeated the VFF 
analysis assuming a LV of radius 6 Mpc both in observa- 
tions and in our 14 simulated LVs. We recovered the same 
circular velocity Vc ~ 35 km s~^ for haloes hosting galaxies 
with AIb < — 12. Thus we are sure that our original sample 
is not influenced by voids which intersect the outer bound 
of the LV. 

However, a signiflcant amount of LV galaxies with mag- 
nitudes Mb ~ —12 and rotational velocity as low as Vrot ^ 
20 km s~^ have been recently observed. A systematically se- 
lected sample of such extremely faint nearby dwarf irregulars 
galaxies h as been observed wi th the Giant Metrewave Radio 
Telescope (|Begum et al.|[200^ ). Dwarf irregular are typically 
dark matter dominated. These objects have surprisingly reg- 
ular kinematics. One of the faintest ones is Camelopar dalis 
B with Mb ~ -10.9 l|Begum, Chengalur fc Hopdl2003h . 

Within the ACDM model mini-voids with sizes Rv > 1 
Mpc cover a volume of about 1600Mpc'' (or about 75 % 
of the sample volume, see Figure [7]). We found a mean 
of 463 ± 85 haloes with 20 < Vc < 35 km s"^ within the 
mini-voids of a simulated LV where the scatter is the la 
fluctuation between the 14 LVs. This number is more than 
an order of magnitude above the observed number of faint 
galaxies in the Local Volume. Thus the ACDM model fails 
to explain the observed local VFF. The predicted overabun- 
dance of low-mass galaxies in low density regions is similar 
to another problem well known in high density regions where 
the ACDM model predicts an order of magnitude overabun- 
dance of low-mass satellites. Typically the more massive 
haloes (30 < Vc < 35km s"^) can be found near the bor- 
ders of the mini- voids whereas the smallest ones tend to be 
in the center. A similar be havior has been found in large 
voids l|Gottlober et al.ll2003l '). 



5.3 RMS peculiar velocities 

Next, the Hubble flow around the LG-like objects is studied 
in an attempt to conflrm that our LV-candidates are indeed 
the numerical counterparts of the observed LV. To this end, 
we calculated the rms radial velocity deviations from the 
Hubble flow, an, both for the observational data and the 
objects in the simulated LVs. As limits of circular velocities 
of haloes used in our analysis, we took the values found 
from the previous void analysis. We used haloes with circular 
velocities Vc > 20 km s~^ in the AWDM and with Vc > 
35 km s~^ in the ACDM sample. But even if we assume for 
the ACDM sample the lower value of Vc > 20 km s~^ the 
result does not change substantially in spite of the signiflcant 

increase in the number of objects. 

Following iTikhonov fc KlvpinI (|2009l ) we calculate an 
in the rest frame of the considered galaxy sample, i.e. we 
correct for the apex motion. We correct the rms velocities 
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Figure 8. The rms radial velocity deviations from the Hubble 
flow an for galaxies with A/g < — 12 in the Local Volume with 
distances from 1 Mpc up to R (red curve with open circles) . The 
estimates are corrected for the apex motion and for distance er- 
rors. Black filled triangles with error bars show the mean an on 
corresponding scales for 14 CDM LV-candidates for haloes with 
Vc > 35km s~^ with \-a scatter. Blue triangles - mean an for 
WDM LV-candidates for haloes with K > 20 km s~^. 

of the galaxies also for distance errors. Depending on radius 
these correctio ns reduce the value of (th by about 5% to 15% 
(cf. Table 3 of IXikhonov fc KlvpinI (120091 )1. Figure H shows 
the observed an in the range from 1 Mpc to R and the 
mean value obtained around the Local Group in the center 
of each of the 14 LVs in the ACDM and AWDM simula- 
tions. In Fig. |8] we compare the an of our numerical LVs 
with the observed one and found a reasonable agreement for 
both models. We error bars show the scatter among the 14 
LVs of the ACDM model. In AWDM we found about the 
same scatter. The excess seen in the observed an at _R = 4 
and 5 Mpc probably reflects the fact that most of the galaxy 
groups of the LV are located within a distance of 3 to 5 Mpc 
from our Local Group. They are responsible for the corre- 
sponding overdensity at these distances. However, we did not 
include this constrain when we searched for LV-candidates 
in simulations for the sake of having better statistics. Not 
surprisingly, the results found from the ACDM and AWDM 
LVs agree pretty well among them and with data. 

Peculiar motions on the scale of a few Mpc reflect the 
mass distribution on these scales and larger ones, hence they 
should not be affected by the small scale damping of the 
primordial perturbations. We conclude that there is no dif- 
ference between the two cosmologies in terms of the internal 
dynam ics of these systems, in good agreement with other 
works IjMartinez-Vaquero et.allboogi ). 



6 FORMING GALAXIES IN HALOES 

As we mentioned above, very faint low-mass galaxies have 
been observed in the local universe with rotational veloci- 
ties below 20 km s~^. As an example, Camelopardalis B with 
Mb ~ —10.9 has a total dynamical mass of « lO*/i~^M0 
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Figure 9. Redshift at which the progenitors of present day 
LV haloes became more massive than3xlO®/i-iM0 as a function 
of present day virial mass (upper axis: number of particles in the 
halo). Filled triangles: ACDM model, open circles: AWDM model. 
Dashed lines: see text. 

measured at the last point of the observed rotation curve 
(|Begum. Chengalur fc Hopdl2003l ). Due to the strong bias 
between dark haloes and galaxies, it is not easy to compare 
these observations with collisionless N-body simulations. A 
proper way of doing it would be by means of full hydrody- 
namical simulations which includes all the relevant physical 
processes that lead to the formation of galaxies inside the 
dark matter potentials. Nevertheless, we can use results from 
such kind of simulations to derive estimates on the formation 
epochs of the galaxies that w ould form in the da rk haloes of 
our simulated local volumes. iHoeft et al.l (|2006l ) have simu- 
lated the formation and evolution of dwarf galaxies in voids 
using high resolution ACDM hydrodynamical simulations 
which include radiative cooling, UV photoionization from a 
cosmic background, star formation from the cooled gas and 
feedback from supernovae. Based on these simulations, they 
derived a characteristic mass scale Mc{z) below which UV 
photoheating suppresses further cooling of gas and therefore 
star formation inside haloes with M < Mc{z). The numer- 
ical results can be well approximated by a exponentially 
decreasing function with redshift in the following way 

M4z) 062 r Ac(0) 1 '''^ 

IQio/i-iM© ^ (l + z)i-23exp (0.08222-1) \Ac{z)j 

(1) 

where Ac{z) represents the redshift dependence of the char- 
acteristic virial overdensity of halos that can be obtained 
from the numerical integration of the equations of the spher- 
ical infall mo del or using the anal ytical approximation ob- 
tained by Bry an fc NormanI (|l998i ) . 

Thus, by comparing the mass of the progenitors of each 
of our LV haloes with Mc{z) we can estimate the typical 
formation redshift for the galaxies in our LVs. Those haloes 
with mass accretion histories that are always above Mc{z) 
have been able to form stars at all times. On the contrary, 
those haloes with mass accretion histories shallower than 
Mc{z) would have not been able to condense baryons in- 
side and they would not host any galaxy. Although these 
results were obtained for ACDM halos, in what follows, we 
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will assume that they can be applied also to our AWDM 
simulations. Since the UV photionization background comes 
from observations and the modeling of the rest of the other 
baryonic processes is the same, the only differences between 
the ACDM and AWDM halos would come from the different 
mass accretion histories. 

The "age" of an object in simulations is often deter- 
mined as the epoch at which half of its mass is gathered. 
Thus massive objects form later. Contrary to this defini- 
tion of formation time, there are others such as the time at 
which progenitors of present day haloes reach a given mass 
threshold. In numerical simulations a natural threshold is 
given by the mass resolution which defines a minimum halo 
mass as the detection limit. Obviously, more massive ob- 
jects reached earlier a given mass threshold and thus, they 
have formed earlier, according to this definition of formation 
time. Here we use the detection limit as an estimate for the 
formation time of the haloes in the simulation and compare 
the evolution of ACDM and AWDM haloes after passing 
this threshold which we take to be the mass corresponding 
to a minimum of 20 particles (i.e 3.2 x 10*/i~^Mq). For the 
cosmological parameters under consideration the character- 
istic mass scale, Mc{z) equals this mass resolution limit for 
haloes in our simulations at redshift z — 3.7 and increases 
up to 6.2 X 10^ h~^MQ at redshift z — 0. Therefore, we are in- 
terested to determine at which redshift Zd the most massive 
progenitor of the redshift zero LV haloes reaches the mass 
threshold (3.2 x lO^/i^^M©) of our simulations. To this end, 
we identify with an FOF analysis all the progenitors of a 
given halo at different redshifts. 

In Figure [9] we show the mean redshift Zd as a function 
of the present day halo's virial mass for all haloes found 
in the 14 LVs of the ACDM and AWDM model. We con- 
sider that for ACDM haloes with more than 100 particles, 
their merging history can easily be followed. In the case of 
AWDM we restrict ourselves to haloes above the limiting 
mass Aliim (see §2). The symbols are centered in the middle 
of the selected mass bins and the error bars show the la 
scatter. The difference between the two models is evident: 
there is a general trend for AWDM haloes of all masses to 
reach this threshold later than the ACDM ones. 

The dashed vertical line in Figure [5] marks the charac- 
teristic mass for redshift zero, Mc(0) = 6.2 X lO^/i"^M0. 
More massive haloes would have form stars all their lifetime 
because their mass accretion histories are always above the 
Mc{z) since za till present. Those haloes with present day 
masses below this characteristic mass would be at present 
practically devoid of gas due to UV photoheating. Within 
the ACDM model the characteristic, Mc, mass increases 
faster with time than the mean halo mass. Therefore, the 
mass of those haloes was at some earlier redshift larger than 
Mc{z) and thus the haloes were able to form stars in the 
past. 

The horizontal dashed line in Figure [5] marks the red- 
shift z* = 3.7 at which M^z*) = 3.2 x lO*/i"^M0 (the 
resolution limit). The redshift of the first identified progen- 
itor of ACDM haloes with Mvir > 2.2 x lO^'/i'^M© is al- 
ways above this line. Therefore, all these haloes have been 
able to cool gas and form stars at least until redshift 3.7. 
Those with masses M > Mc(0) will have formed stars up 
to present. The smaller haloes, with M < Afc(O) would 
have formed stars only until Mc{z) became larger than the 
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Figure 10. The differential number fraction -jY~"4^ of progen- 
itors of haloes which became more massive than 3 X 10^/i~^Mq 
at redshift z^- Filled circles: haloes with Vc < 35 km s~^ at z = 
inside mini-voids, open circles: haloes with Vc < 35 km s~^ out- 
side mini- voids, triangles: haloes with Vc > 35 km s~^. Error bars 
are la deviations. 



progenitors halo mass (somewhere between Zd and present). 
Haloes with Mvir < 2.2 x W'^'h~^MQ lost most of their gas 
before z* = 3.7. Therefore, they can host only an old popu- 
lation of stars. 

In sharp contrast with the ACDM case, the AWDM 
haloes with Mvir < Mc(0) are always below the horizontal 
dotted line. Thus their mass was always smaller than Mc{z) 
and they were unable to condense baryons and form galax- 
ies at all. This points towards a potential problem of the 
AWDM model, namely an insufficient formation of dwarf 
galaxies. Since Mc{z) drops down very fast for 2 > 4, in 
principle they could have been able to cool gas and form 
stars at earlier times, resulting in a population of faint red 
dwarfs. 

Since haloes form much later in the AWDM scenario 
they can lose their gas much easier due to the photoheating 
of the cosmic ionizing background than the ACDM haloes 
of the same present day mass. Therefore, the formation of a 
sufficient number of low mass galaxies could be a problem in 
this model. To study this in full detail, one needs to perform 
AWDM simulations with much higher mass resolution than 
the one presented here. 

In Figure [TO] we compare the formation history of 
ACDM haloes with circular velocities Vc < 35 km s~^ lo- 
cated inside the mini- voids (filled circles) with haloes in the 
same range of circular velocities, but located outside the 
voids (open circles). We show the differential number frac- 
tion of progenitors of haloes which became more 
massive than 3.2 x 10^h~^MQ at redshift Zd which has the 
same meaning as in in Fig. (9] We used haloes with masses 
M > 2 X W^h'^MQ. Therefore, the mean Zd of Fig.lTOlcor- 
responds to the second filled triangle of Fig. O As can be 
seen in the figure, regardless of environment, these haloes 
have very similar formation histories. For comparison, we 
also show the results for the more massive haloes in the LVs 
(triangles: Vc > 35 km s~^) which obviously passed much 
earlier the threshold Zd- The fact that we do not see strong 
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differences in the formation epochs of dwarf haloes located 
in structures and those sitting in voids may have profound 
consequences for the ACDM model. 

One of the proposed scenarios to solve the overabun- 
dance problem assumes that UV-photoionization evaporates 
gas from dwarf haloes in voids and therefore quenches star 
formation and leaves voids empty. If, however, dwarf haloes 
in voids and the surrounding structures form roughly simul- 
taneously then photoionization should have suppressed the 
formation of all these dwarf galaxies regardless of environ- 
ment. As shown above the ACDM model predicts an order 
of magnitude more low-mass haloes than observed. If pho- 
toionization suppresses the formation of galaxies in more 
than 90 % of the small haloes, it is difficult to explain why 
none of the remaining is found in the voids. Therefore, the 
observed small number of dwarfs represents a similar prob- 
lem for the ACDM model as the missing satellites. 

However, the observed dwarfs could also be a poten- 
tial problem for the AWDM model. According to the late 
formation time of the small mass haloes they spend almost 
all their lifetime having masses well below the character- 
istic m ass Mc{z). Following the arguments of iHoeft et al „ 
l|2006l ) it would be difficult to explain how they could have 
formed stars at all. However, we should also note that the 
charact eristic mass has been derived based on ACDM simu - 
lations l|Hoeft et al.ll2006l : lOkamoto. Gao. fc Theunsll2008h . 
To get a final answer for halos in the WDM scenario, proper 
gas-dynamical simulations including star formation should 
be performed. Such simulations should have extremely high 
resolution to address also the problem of fake halo forma- 
tion. 



7 DISCUSSION AND CONCLUSIONS 

We u sed an updated version of the iKarachentsev et al.l 
l|2004l ) sample of galaxies to compare the distribution and 
motions of galaxies in the Local Volume with that of dark 
matter haloes in ACDM and AWDM simulations. About 
40 mini- voids with sizes in the range from 0.5 Mpc to 
4Mpc have been detecte d in the observational sample. 
iTikhonov fc KlvpinI ()2009l ) have shown that the spectrum 
of void sizes is stable with respect to variations of the sam- 
ple and to uncertainties of distances to individual galaxies. 

From our ACDM and AWDM simulations we have se- 
lected t he Local Volume candidat es following the criteria 
used bv lTikhonov fc KlvpinI l|2009l) . In many respects these 
LV-candidates look very similar to the observed Local Vol- 
ume: they have comparable number of bright (massive) ob- 
jects and a similar density contrasts on 8Mpc scale. More- 
over, most of the LV-candidates show the same typical inter- 
secting filaments sometimes nearly aligned in planes which 
are quite similar to what one observes in the Local Volume 
- intersecting filaments form a structure which is part of 
the Local Supercluster and goes through the Local Group 
(sometimes it is called "Local Sheet"). The 14 Local Volume 
candidates in ACDM and AWDM simulations show only 
little scatter in the volume fraction occupied by mini- voids. 
Thus it is well suited to discriminate between different mod- 
els. To test further our LV-candidates we have calculated 
the rms radial velocity deviations from the Hubble flow and 
found no difference in the dynamics of the haloes. Both in 



the ACDM and AWDM models an is in agreement with the 
observations. 

In Sect. 13.11 we argue in favor of an association of the 
haloes maximum circular velocity to the observed rotational 
velocity. Our argument is based on a comparison of the ro- 
tational velocities observed in a certain distance from the 
center and the circular velocities of dark matter haloes in 
the same distance. The velocity profiles of the dark matter 
haloes in are essentially flat at the distance where the rota- 
tional velocities of dwarfs are measured. Thus the errors are 
small if one associates the measured rotational velocity to 
the peak circular velocity. The exact relation between both 
depends on the gasdynamical processes in the galaxy which 
we cannot consider in the large volume under discussion. 
One might argue that the peak velocity increases due to adi- 
abatic contraction so that the WDM model would be even 
more favored. However, these dwarf haloes are expected to 
be dark matter dominated and adiabatic contraction is ex- 
pected to be less important. One might argue that by some 
processes the rotational velocity in the dwarf haloes becomes 
smaller but we cannot see a good reason to associate Vrot 
with significantly more massive haloes which would be nec- 
essary to fit the observational data. 

Within the ACDM model the cumulative volume frac- 
tion of the LV occupied by mini-voids matches observations 
only if one assumes that objects with Vc > 35 km s~^ define 
the local mini- voids and these objects host ga laxies brighter 
than Mb = -12 l|Tikhonov fc KlvpinI [20091 ) . The ACDM 
model predicts almost 500 haloes with 20 < Vc < 35 km s~^ 
within the mini-voids in a LV. At present 10 quite isolated 
dwarf galaxies have been observed in the Local Volume with 
magnitudes —11.8 > Mb > —13.3 and rotat ional velocities 
Vrot < 35 km s"^ (|Tikhonov fc KlvpinI [2OO9I '). This number 
is more than an order of magnitude smaller than the pre- 
dicted number of low-mass haloes and thus it points to a 
similar discrepancy as the well known predicted overabun- 
dance of satellites in Milky Way sized ACDM haloes. The 
reason can be understood from the ACDM mass function 
at the low-mass end which is much steeper than the ob- 
served luminosity function at the low luminosity end. This 
causes both the overabundance of satellites and the fact, 
that the observed mini-void distribution can be explained 
only if haloes with circular velocities larger than 35 km s"'^ 
host galaxies. In case of satellites one could explain the mis- 
match between observations and predictions by processes 
like ram-pressure and tidal stripping which have not been 
taken correctly into account in simulations. On the contrary, 
the dwarfs in the Local Volume are all situated in the same 
environment of relatively low density. Therefore, the exis- 
tence of only a small number of such low luminosity and 
low mass dwarf galaxies in the LV is even more difficult to 
explain than the missing satellites. 

In case of the AWDM model the observed spectrum of 
mini- voids can be explained naturally if dark matter haloes 
with circular velocities larger than ~ 15 — 20 km s~^ host 
galaxies. The observed small number of slowly rotating iso- 
lated galaxies is not in contradiction with this assumption 
since a small number of low-mass haloes are also formed 
in the AWDM model. Thus the AWDM model solves the 
problem of overabundance and the spectrum of mini-voids 
in a distribution of low-mass haloes is in agreement with 
observations. However, we have also shown that due to the 
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late formation of haloes in the AWDM model a substan- 
tial part, if not all, of the low-mass haloes could have failed 
to form stars. Thus they would have not been able to host 
any visible galaxy, as it is necessary for the model to agree 
with the sizes of mini-voids in the Local Universe. Never- 
theless, to study their distribution and properties in detail 
one needs simulations with much higher resolution so that 
one can distinguish between the real low-mass haloes and 
the numerical artifacts. Moreover, full radiative hydrody- 
namical simulations with star formation are also needed in 
order to properly address the galaxy formation process in 
these halos. Overall we confirm earlier claims that concern- 
ing the formation of small scale structure the AWDM model 
can be considered as a promising alternative to the ACDM 
model. However, more detailed investigations are necessary 
to study to which extend star formation in low-mass haloes 
is suppressed in the AWDM model. 
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